Regular atomic narrowing of Ni, Fe and V nanowires resolved by 2D correlation analysis 
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We present a novel statistical method for the study of stable atomic configurations in breaking nanowires 
based on the 2D cross-correlation analysis of conductance versus electrode separation traces. Applying this 
method, we can clearly resolve the typical evolutions of the conductance staircase in some transition metal 
nanojunctions (Ni, Fe, V) up to high conductance values. In these metals our analysis demonstrates a very well 
ordered atomic narrowing of the nanowire, indicating a very regular, stepwise decrease of the number of atoms 
in the minimal cross section of the junction, in contrast to the majority of the metals. All these features are 
hidden in traditional conductance histograms. 

PACS numbers: 73.63.Rt,81.07.Lk,73.23.-b,61.46.Km,07.05.Kf 
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The study of electronic transport on the single atom or 
single-molecule scale have revealed several interesting phe- 
nomena due to the interplay of the atomic granularity of mat- 
ter and the quantum nature of conductance l li]. The conduc- 
tance of stable, frequently occurring atomic configurations 
or molecular contacts formed during the rupture of metallic 
nanowires can be identified using peaks in the conductance 
histograms constructed from a large number of conductance 
versus electrode separation traces [1-7]. Individual conduc- 
tance traces exhibit several plateaus up to tens of the conduc- 
tance quantum unit. Go = 2e^//i in accordance with simula- 
tions predicting a stepwise decrease of the number of atoms 
in the minimal cross-section (MCS) of the junction (jsl 
In the most well studied system, gold the discreteness of the 
MCS can be reflected by numerous peaks in the histogram 
lUt]. However, in the majority of the cases the histograms only 
show a single or a few peaks fH; thus the knowledge about 
the details of contact rupture is very limited. This is especially 
true for transition metal nanojunctions, where the conduction 
through low-symmetric d orbitals introduces an extreme sen- 
sitivity to the contact geometry giving rise to a broad variety 
of conductances for contacts with the same MCS lIsilToll. 

For a deeper understanding of precise atomic arrangements 
in the junction the conductance histograms must be supple- 
mented with further experimental input, as well as by detailed 
theoretical ab-initio simulations. It was already recognized 
that additional statistical analysis of the conductance traces 
can provide fundamental information, like the discovery of 
monoatomic chains by plateaus' length analysis flT], or the 
deeper understanding of single-molecule junctions by recently 
developed statistical methods 

m. 

In this Letter we introduce a novel characterization method 
based on the 2D cross-correlation histogram (2DCH) analy- 
sis of conductance traces lll2ll in some analogy to 2D corre- 
lation spectroscopy (COSY) in magnetic resonances [l?]. In 
the field of mesoscopic physics it is already known that cor- 
relation analysis provides significantly new information com- 
pared to the measurement of mean values, which was demon- 



strated by the success of quantum noise studies lfl4ll . Imple- 
menting these ideas we demonstrate that the statistical cross- 
coiTelation analysis of conductance traces opens a possibil- 
ity to answer several questions concerning the atomic-scale 
nanocontact formation dynamics. Applying the correlation 
technique we study the contact evolution in some transition 
metal nanojunctions (Ni, Fe, V), demonstrating that the con- 
ductances of plateaus related to different atomic configura- 
tions show correlated fluctuations. In the conductance his- 
togram these fluctuations introduce a strong smearing; how- 
ever, the 2DCH clearly resolves the alternative evolutions of 
contact conductance demonstrating a very regular stepwise 
change of the MCS starting from rather high conductance 
values. In comparison, this very regular naiTowing of the 
nanowire is not detected in the majority of the metals. 

2D correlation histograms. — We begin the discussions 
using our result for Au nanowires chosen as a well charac- 
terized reference system. The measurements were performed 
on high purity samples applying the mechanically control- 
lable break junction technique at liquid helium temperature 
0151]. Typical conductance vs electrode separation traces fea- 
ture some similarities in the general character and differences 
in the details due to the stochastic nature of contact rupture 
[Fig. m^)]- A conductance histogram for a chosen trace, r 
can be constructed by counting the number of data points in 
discrete bins, i along the conductance axis: Ni{r) [Fig.[Tlb)]. 
The histogram for the whole data set including several thou- 
sand traces [Fig. [TJc)] is defined by the average number of 
data points in the different bins: Hi — {Ni{r))r- 

Though conductance histograms have a fundamental role in 
the study of atomic and molecular nanojunctions, it is evident 
that histograms grab only partial knowledge from the mea- 
sured data. Significantly more information can be collected 
by studying the cross-correlations between different atomic 
configurations. The cross-correlation function between bin i 
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FIG. 1: (Color online) Measurements on Au junctions: (a) represen- 
tative conductance vs electrode separation traces; (b) histogram for 
the middle trace on panel (a); (c) histogram for the whole data set in- 
cluding 5000 traces; (d) 2D correlation histogram, dj for the same 
data set. (e) Illustration for slant plateaus and neighboring jumps 
respectively being responsible for positively correlated areas at the 
diagonal of the correlation plot and the neighboring negatively cor- 
related region. 



and j can be defined as: 
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where SNi{r) = Ni{r) — {Ni{T))r- This function gives zero 
if Ni{r) and Nj{r) are independent from each other. If the 
occurrence of two plateaus near bins i and j is correlated - 
i.e. either both of them or none of them appear - the correla- 
tion function takes positive values. Contrary, it takes negative 
values if two plateaus are anticorrelated; i.e., the occurrence 
of a plateau near bin i excludes the occurrence of another one 
at bin j. By definition the correlation is limited to the range 
— 1 < Cij < 1, it is symmetric, Cij = Cj^i, and the diagonal 
always gives Ci^i — 1. The cross-correlations can be visu- 
alized by a 2D color plot, where the two axis correspond to 
bins i and j and the color-scale shows C,;.j with warm colors 
(yellow-red) for positive coiTelations, cold colors (blue-black) 
for negative coiTelations, and the areas where the coiTelation 
is zero within the accuracy of the method is marked by green. 

In Fig.lHd) we present the 2DCH for Au nanoj unctions. An 
evident feature observed on the 2DCH is a finite width posi- 
tively correlated region near the diagonal, and a wider nega- 
tively correlated region near both sides of the diagonal. This 
feature has a simple reason: if on a certain trace a plateau 
is observed near bin i (i.e. the number of points in bin i is 
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FIG. 2: (Color online) Measurements on Ni junctions, (a): Bottom 
area graph shows the conductance histogram for the whole data set 
including 5000 traces. The line graphs show conditional histograms 
for different subsets of the traces (see text). Horizontal dashed lines 
show the zero values of the histograms, vertical dashed lines indi- 
cate the peak positions in the middle conditional histogram, (b) 2D 
correlation histogram for the whole data set. (c) Representative con- 
ductance traces, (d) Illustration for the correlated shifting of plateaus. 



larger than in average) then, assuming a plateau with finite 
slope, also the neighboring bins will have larger number of 
points than in average, introducing positive correlation. Fur- 
ther away from bin i jumps are observed on the trace, intro- 
ducing negative correlation [see demonstration in Fig. [TJe)]. 
For even larger conductance differences no significant corre- 
lations are detected. 

Correlation analysis ofNi, Fe and V junctions. — We have 
performed 2DCH analysis of some transition metal nanojunc- 
tions: Ni, Fe, V. Here, we present and discuss our results on Ni 
junctions in detail [Fig.|2l, while the experiments on Fe and 
V exhibiting similar features are demonstrated in the supple- 
mentary material flE], and discussed in detail elsewhere [1^. 
The typical histogram of Ni, Fig. |2a) exhibits two peaks at 
G = 1.7 Go and G = 3 Go, but at higher conductances no 
clear features are resolved, in accordance with previous low- 
temperature studies [i TtI [Tsll . In contrast, individual conduc- 
tance traces [Fig. Efc)] show clear plateaus up to tens of the 
conductance quantum unit. The conductance of the plateaus 
exhibit considerable variation from trace to trace which causes 
a complete smearing of the histogram at higher conductance 
values, and only two low-conductance peaks survive. 

We present the 2DCH for Ni contacts in Fig. |2|b), where 
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FIG. 3: (Color online) Simulations for stripes, (a) Experimental data 
on Ni with two simplified fits, (b) More advanced simulations. (See 
text and supplementary material for details). 

similarly to Au junctions, an anticorrelated region is seen be- 
side the diagonal, however further away several strong posi- 
tively correlated stripes are observed, which are almost paral- 
lel with the diagonal. The regions positively correlated at the 
coordinate correspond to the simultaneous appearance 
of plateaus both at bins i and j; thus the first stripe close to 
diagonal defines the probable conductance pairs of neighbor 
plateaus. The presence of stripes indicates that the conduc- 
tances of plateaus change in a correlated way: if the conduc- 
tance of a certain plateau shifts upwards or downwards from 
trace to trace due to the changes in the contact environment, 
the conductances of neighbor plateaus are shifted in corre- 
lated way leaving the differences in conductance similar [see 
demonstrative cartoon in Fig.|2jd)]. 

To better resolve the nature of the correlation, in Fig. |2a) 
we plot separate histograms for selected traces that have a 
plateau at a certain conductance - the top, middle and bot- 
tom curves correspond to conductance traces which have the 
final plateau near G — 1.6, 1.7, 1.8 Go, respectively. Pre- 
cisely, these so called conditional histograms are constructed 
from those traces, r', for which the number of data points in 
a certain bin, j is larger than its average value: iJ^"""* — 
{Ni{r')\SNj{r') > 0)r'. With this objective selection crite- 
rion the traces having a long enough plateau near bin j can be 
sorted out with reasonable certainty. These conditional his- 
tograms exhibit up to several peaks, which is the most pro- 
nounced for the middle graph, where the final plateau before 
complete rupture is fixed at the position of the first peak in the 
total histogram. For traces with somewhat smaller or larger 
conductance of the last plateau, all the peaks are respectively 
shifted to lower or higher values. This demonstrates that the 
nearly featureless total histogram builds up from partial his- 
tograms with several well-defined peaks. 

On the magnetic samples, Ni & Fe control measurements 
were performed in high magnetic field. A magnetic origin 
of the reported features can be excluded by observing similar 
conductance histograms and 2DCHs at B = 12 T. 

Correlated shifting of plateaus. — Both the 2DCH and the 
conditional histograms for Ni nanojunctions clearly show the 
correlated shifting of plateaus from trace to trace. For a better 
understanding of this phenomenon we have performed simple 
simulations for the observed stripes, which are demonstrated 



in Fig. |3] and discussed in more detail in the supplementary 
material [15]. The resistances of plateaus (R — G^^) can be 
modeled as the resistance of some base configurations being 
modified due to the random variations in the contact neigh- 
borhood, Rff = i?^^'" + ARn^r, where n and r label the 
different plateaus and traces, respectively. In Fig. [Sfa) we 
present the results of two simple models, in which the vari- 
ation of the resistance is random for the different traces, but 
on a certain trace it is either constant for all the configurations 
{ARn^r = ARr, curvcs above diagonal), or it is proportional 
to the base resistance (Ai?„.,. = • R^'^'', curves below di- 
agonal). Fig. Ob) shows that finite width stripes very much 
resembling the experimental data can be obtained by more de- 
tailed simulations generalizing the latter model llSll . 

Various physical mechanisms can influence the conduc- 
tance along several steps of the evolution of the junction, and 
thus lead to a correlated shifting of the plateaus. These can be 
the trace to trace variations in the crystallographic orientations 
in the vicinity of the contact, changes in the contact geome- 
tries, or scattering on extended defects [19] being pinned to 
the contact region. 

Discreteness of the minimal cross section. — Our results 
show some analogy to previous studies on work hardened Au 
and Al samples [2, 3], where numerous peaks were observed 
in the conductance histogram reflecting the quantized number 
of atoms in the MCS of the junction. In Au, where symmet- 
ric s orbitals contribute to the conduction, the conductance is 
mainly proportional to the MCS area [2], which explains that 
due to the discreteness of the MCS the conductance of the 
junction also prefers some discrete values. In Al the observa- 
tion of similar features is already somewhat surprising due to 
the role of p orbitals in the conduction [i3J . In transition met- 
als already up to six conductance channels can contribute to 
the conductance through a single atom [1, 8, 10], from which 
the conduction through low-symmetric d orbitals is extremely 
sensitive to the directions of the bonds between the atoms. 
Therefore, contacts with the same MCS but different geome- 
try (e.g. different crystallographic orientations) have a broad 
variance of conductances [8] leading to a strong smearing of 
the histogram. Our measurements on Ni, Fe and V junctions 
show that - in spite of the large number of plateaus in the con- 
ductance traces - the histogram is indeed smeared, and the 
peaks are absent at higher conductance values. However, the 
2DCHs can resolve the correlated shifting of plateaus giving 
a unique opportunity to follow the statistically relevant evo- 
lutions of the conductance staircase: once we fix the starting 
conductance of a high conductance plateau, it very well de- 
termines the subsequent plateau conductances during the rup- 
ture, as demonstrated by the arrows in Fig. I3b). These well- 
defined sequences of plateau conductances, and the numerous 
peaks in the conditional histograms demonstrate a very reg- 
ular, well ordered rupture, indicating a well-defined discrete 
sequence of atomic configurations during the narrowing of the 
nanowire. 

In contrast to the equidistant peak positions in the his- 
tograms of work hardened Au [2], the distances between the 
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peaks in the conditional histograms for Ni increase with in- 
creasing conductance. In the 2DCH the average difference in 
conductance between neighboring plateaus - defined by the 
horizontal (or vertical) distance of the first stripe from the di- 
agonal [Fig.|2|b)] - accordingly shows a smooth increase. The 
low G value of the conductance jump, AG « 1.6 — 2 Go is 
fully consistent with the removal of a single atom from the 
MCS in each step. At higher conductances AG increases up 
to 3.5 Go, which can still be consistent with removing a sin- 
gle atom with several open conductance channels, assuming 
that the average conductance per atom is increasing with in- 
creasing junction size. Indeed, for d metals it may happen 
that low coordinated atoms on the circumference of the junc- 
tion give smaller contribution to the conductance than fully 
coordinated atoms in the contact center due to the role of d or- 
bitals oriented with large angle to the junction axis. Based on 
these simplified considerations the atom by atom narrowing 
of the junction is a possible explanation for the experimen- 
tal observations. A regular, locally crystalline packing of the 
electrodes with slips of the lattice planes by atomic distances 
could also account for the observed regular rupture. However, 
an unambiguous microscopic explanation of the observed reg- 
ular narrowing can only be given in comparison with detailed 
density functional theory calculations. 

Contrary to Ni junctions our measurements on Au did not 
show any stripes in the 2DCH [Fig. [lid)]. There is a sig- 
nificant difference in the mechanical behavior of two met- 
als, which is clearly demonstrated by the traces presented in 
Figsina) and He). In Ni junctions the conductance traces 
are very well ordered with well-defined conductance jumps 
[Figl2c)], whereas in Au - though several plateaus are ob- 
served - the conductance jumps shows a rather broad vari- 
ety of amplitudes [FiglTJa)]. This demonstrates, that - al- 
though the numerous peaks in the histogram of Au reflect the 
discreteness of the minimal cross section - the conductance 
histogram cannot tell how the MCS changes during a single 
conductance step. According to our measurements, in Au the 
changes of the atomic configurations are not regular, very fre- 
quently either multiple atoms are removed from the MCS, or 
just smaller rearrangements appear with AG < 1 Go, which 
smears the structures in the 2DCH. In Ni, however, the am- 
plitude of the conductance jump is a very well-defined func- 
tion of the junction size, as shown by the robust stripes in the 
2DCH. This points to a very well ordered mechanical evolu- 
tion of the junction, and indicates that the junction follows a 
well-defined sequence of atomic configurations. To check the 
relevancy of this unique behavior, we have performed mea- 
surements on a broad variety of metals flBIl . We have found 



that this regular atomic narrowing is absent in the majority 
of the metals, including the noble metals Au, Ag, Cu; the p 
metals Al, Sn, Pb; and transition metals Ta, Nb. In Pt and Pd 
junctions we have detected some signs of regular narrowing, 
but it was less pronounced than in Ni, Fe and V. 

Conclusions. — We have demonstrated, that several fea- 
tures of atomic-scale contact formation dynamics - being hid- 
den in traditional conductance histograms - can be resolved 
by the 2D cross-correlation analysis of the conductance traces. 
Applying this method we were able to identify a very reg- 
ular atomic narrowing during the rupture of Ni, Fe and V 
nanowires, which is absent in the majority of the metals. 
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